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The universal power law tails of single and multi-particle time correlation functions are derived 
from a unifying point of view, solely using the hydrodynamic modes of the system. The theory 
applies to general correlation functions, and to systems more general than classical fluids. Moreover 
it is argued that the collisional transfer part of the stress-stress correlation function in dense classical 
fluids has the same long time tail ~ £ _1_d / 2 as the velocity autocorrelation function in Lorentz gases. 

PACS numbers: 05.20Dd kinetic theory 

05. 40. -a Fluctuation phenomena, random processes, noise, and Brownian motion 

The long time tail (LTT) of the velocity autocorrelation function (v x (t)v x (Q)) ~ t~ d / 2 in d— dimensional 
fluids in thermal equilibrium is the prototypical example of time correlation functions in fluids showing 
universal power law behavior, independent of the details of the inter-particle interactions The goal of 

this paper is to present, from a unified point of view, a quantitative description of the asymptotic decay of 
single and multi-particle time correlation functions (TCF's). The starting point is the mode coupling theory 
for classical fluids in 2], which is extended to more general systems, and more general TCF's. 

What about the asymptotic decay of other single particle properties, that are coupled to a conservation 
law? Consider the single site (n = 0) correlation function (vo(i)vo(O)), initialized in a state of thermal 
equilibrium, which state may or may not be maintained by the dynamic evolution of the system. The 
system considered has fixed spins or velocity vectors v n at the sites n = {m, n 2 , n^} of a d— dimensional 
cubic lattice with nearest neighbor (n.n.) interactions conserving the total magnetization or momentum, 
^ n v n (i). For convenience we take periodic boundary conditions. Such models are Glauber's model with 
Ising spins at zero temperature 0, or Lattice Gas Cellular Automata with b velocities per site associated 
with the links , or granular fluid models with a discrete or continuous velocity vector assigned to each site 
0, 0? ] . In the following we consider as an example an exactly soluble model of a granular fluid on a lattice 
whose LTT's are also covered by the results, to be derived in this letter. 

Here a pair of n.n. sites interacts at a rate kq such that v n and v m are each replaced by their arithmetic 
mean. Then the equation of motion for the mean value, (v n ) = V n , can be written in appropriate units of 
length and time as a discrete diffusion equation, i.e. 

dV n /dT = i Ea[ V n+a - V n ] = ±AV n (1) 

Here r = 2Dt/a 2 is the rescaled time, D = n a 2 the diffusion coefficient, a the lattice distance, a runs 
over n.n. sites, and A is the Laplace operator on a d— dimensional cubic lattice. The exact solution 0,0] 
of this equation is V n (r) = J2 m V n -m(0) T\'i=i[Irn i (T)e~ T }, where I m (r) is the modified Bessel function. 
The space-time correlation function (v nx (t)vo x (0)) satisfies the same equation as (v nx ), where (• • • ) is an 
ensemble average over an arbitrary initial state. In the sequel all TCF's refer to initial states in thermal 
equilibrium. They are not necessarily spatially uniform. The single site TCF, Co(i) = (vo x (t)vo x (0)) / (vq x ) , 
follows by setting n = 0, and yields the exact solution, 

Co(t) = (e- T Io(r)) d ~[2KT]- d ' 2 (2) 

where the asymptotic equality gives its LTT. Here Co(t) is essentially the return probability of initial mo- 
mentum to its points of origin. The same exact results apply to the spin TCF in the Glauber model at zero 
temperature, and to the energy auto-correlation function C e (t) — (Se n (t)Se n (o)) / ((Se) 2 ) in case the total 
energy ^\ Ei(t) is conserved, where Se n — e n — (e). This function is again a return probability. 

In all previous models particles have only n.n. interactions, and are placed on lattice sites. These restric- 
tions will be removed. When the particles are moving, as in fluids, one may also consider the tagged particle 
fluctuations, j% = v 2 — (v 2 ), or more generally j 2 k — v 2k — (v 2k ), and j2fc+i = v 2k v x , ... with k = 0, 1, 
and even jo = 5e in case the tagged particle has an additional internal degree of freedom with a discrete or 
continuous energy e. Here the total energy of the system is conserved, and we are interested in the LTT of 
these TCF's. Such systems are classical fluids DPD fluids (Dissipative Particle Dynamics), which are 
mesoscopic models 0, with particle, momentum and possibly energy conservation, or the LBE (Lattice 
Boltzmann equation) method which lacks energy conservation. The DPD fluid and LBE method can 
be considered as pre-averaged versions of respectively the Liouville equation and the dynamic equations of 
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In the DPD fluids the evolution equations for {ej(t), Vj(t)} are formulated as coupled Langevin equations 
with dissipative and stochastic pair- interaction terms, having a finite interaction range r c , and satisfying the 
fluctuation dissipation theorem. The stochastic interactions contain in general multiplicative white noise 10] . 
One may also quench the translational degrees of freedom, and freeze the particles at fixed random positions, 
or in fixed periodic lattice configurations. In the former case a final average over quenched configurations 
has to be performed. Here we consider only the quenched DPD fluid with internal energy states, referred 
to as DPD solid |l0l where dissipative and random interactions conserve the total internal energy. 
This quenched model of " interacting heat particles" is in several respects the dual model of the overlapping 
Lorentz gasflfjj. 

The purpose of this paper is to study the LTT's of single- and multi-particle TFC's, ( Jb(t) Jft(O)), in 
different models, and away from percolation and other critical points, using mode coupling theory. We start 
with fluid models. Here the LTT's can be concisely summarized in the Kadanoff-Swift formula [2j], 

(Mt)j b (o)) * I4EaJ^1 2 ^ + < )4 (3) 

where J k = V r_1 ^ k . In the thermodynamic limit J k is replaced by (27r) _d J dk. The A/x-sum extends 
over the complete set of hydrodynamic modes {a k } of the system. The amplitude A b ^ = (</&|ak a -k) = 
V~ x ( Jfca^O-k) represents the component of J b parallel to a product of two modes. These modes are linear 
combinations of Fourier transforms of conserved densities, and are normalized as (flkl a k) = ^V- The 
scalar products are defined as {A^B^) — V -1 (A^B^), averaged over an equilibrium ensemble. All TCF's 
considered here approach zero for large times, which implies that (Jb\a£) = 0. 

In a classical fluid there are (d— 1) diffusive shear modes (A — rji : i = 1,2, • • • , d—1), being the components 
of the tranversal flow field Uk± = Uk kk ■ Uk, with decay rates z^ = —uk 2 and shear viscosity r\ = pv, and 
one diffusive heat mode with decay rate = —Dxk 2 and heat diffusivity Dt- In addition there are 
two damped propagating sound modes (A = a = ±) with decay rate z% = —icrcok — ^T s k 2 , where Co is the 
speed of sound and T s the sound damping constant. For the explicit expressions of fluid modes a k in terms 
of conserved densities, and the calculation of the dominant LTT's we refer to Q. 

A special case are the single particle TCF's (jb(t)jb(0)) in fluids. Here the mode coupling formula nec- 
essarily involves the Fourier mode of the tagged particle density, i.e. the self diffusion mode, a k = n k = 
exp[— ik • ri], where i = 1 is the label of the tagged particle. Consequently the following replacements have 
to be made in ©: /i — > s and ^X)a^ ~~ > Sa- Moreover, if one of the labels equals s, then the amplitude 
Ab S = C/fcl a k a -k) = 0fc a k a -k) without a factor l/V. To extract the long time behavior from © we change 
the integration variable k to q/v^, and take the long time limit, using a k — > Oq and n k — > — 1. So, the 
amplitude of the single particle function simplifies to = (j/,|ao), and depends only on the direction k of 
the wave vector k. The final result for the LTT becomes, 



(j b (t)j b (0)) * £ [{j b H)] 2 /[MDx + D)t\ d ' 2 (4) 

A 

where only diffusive modes (A ~ (r)i,H)) contribute. The overline indicates an average over the solid 
angle k. As the calculations of A^ s are similar to those in [2b], we simply quote the final results for the 
LTT in the equilibrium TCF's of classical fluids, C b (t) = (jb(t)jb(O)) / with b = {1,2,3,4,...}, and 

jb = (v x ,v 2 - (v 2 ),v 2 v x ,v 4 - (v 4 )), i.e. 

Ci(i) ~ [{d-l)/nd][4ir{is + D)t}- d / 2 

C 2 (t) ~ (C v /nC p )[An(D T + D)t]-^ 2 (5) 

with C p the specific heat at constant pressure, and C° = \dks- One similarly shows that odd-in-w, c.q. 
even-in-?; correlations are at large times proportional to C\(t) and (^(t), i.e. Cz(t) ~ \{d + 2)/(d + 4)}C\{t), 
d{t) ~ [(rf + 2)/(d + 3)lC , 2(i), etc for k = 5, 6, • • • . The result for the VACF in hard sphere or Lennard-Jones 
fluids is well known [l|,|2j. The remaining ones are new. The above results also disprove the misconception 
that momentum conservation is necessary for the existence of LTT's in classical fluids. The results in (QJ-lJSJ) 
for classical fluids also apply to DPD fluids with energy and momentum conservation 0. In DPD fluids 
without energy conservation only the relations for odd- A: values apply because the heat mode is absent, 
and the system is thermostatted instantaneously. 

Next we want to extend these results to the random DPD solid. The mode coupling results and 
(0} for TCF's in fluids do in general not apply to systems with quenched disorder jl^]. There they apply 
only to TCF's C(t\X), calculated in a non-uniform equilibrium state, corresponding to a single quenched 
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correlation, C(t) = (C(t\X)), a subsequent average over all quenched configurations has to be performed. 
In general the X— dependence of A X)i (X) and Dt(X) is not known explicitly. A more phenomenological 
derivation of a mode coupling formula for diffusive systems with quenched disorder has been presented in 
Ref.^3|j an d one may try to extend that method to the heat conducting random DPD solid. However, we 
will not follow that route here, but investigate only those special cases for which the LTT of C(t\X) can be 
determined explicitly from J2J. 

Consider first the single particle energy correlation C e (t\X) = (<5ej(i)<5ej(0)) x m the DPD solid, where the 
heat mode, a^(t) = a£f (0) exp[— ifc 2 D<r(X)], is the only slow macroscopic mode, and = e^/ J (ek|et). 
In the long time limit (where k — > 0) the relevant amplitude A Hs (X) = (5eia^af_ k ) x — (Seia§)x = 

\/V((8e) 2 ) /N where a k = exp[— ik • ri] is a frozen mode with z| = 0. Inserting these results in Q, and 
performing the X— average yields the LTT, 

v ( i \ d/2 1 

Ge{ - t) ~ N \^Drt ) = p[4nt*] d / 2 (6) 

Here (Dt(X)) = Dt, and fluctuation corrections of relative order ((5Dt(X)) 2 ) / Dt have been neglected. 
Moreover t* = Dt/r 2 and p = Nr 2 /V is proportional to the mean number of particles inside an interaction 
sphere of radius r c . If the particles are put on a lattice (no disorder), the LTT in C e (t) applies as well to 
the lattice version of the DPD solid. Both on-lattice and off-lattice computer simulations of the LTT in C e 
are in excellent agreement ^lj with the theoretical prediction © . 

Next we study the N— particle correlation Cqit) = (Q x (t)Qx(0)) /V . The random DPD solid sustains a 
microscopic heat flux, Q x = Qd + Qr, with a dissipative (D) and a random part (R), as given explicitly 
in ^J. We start with Qd{1) = Ko 12i<j w ( r ij •) r ij,x( e j{t) — Ej(t)). This sum of dissipative pair interactions 
represents the instantaneous exchange of energy over a distance through the interactions (collisional 
transfer). The model parameter kq represents the interaction frequency, and the range function is a step 
function, vanishing for r > r c , and normalized such that J drw(r) — rf. Consequently Ko( £ j — e i) w ( r ij) 
is the rate of energy transfer between the interacting pair (ij) with r»j < r c . The corresponding TCF is 
Cz>(i) = V^ 1 (Q o(t)Q d{0)) ■ Next we consider the stochastic part of the heat flux Qr, and the corresponding 
correlation Cn(t) — V~ x (Q_r(^)Q_r(0)) and cross-correlations CDn{t) and Cjn)(t). For our present purpose 
it is sufficient to note that Qr is linear in the Langevin force, Fij(t) = —Fji(t), with coefficients depending 
on e,; and ej (multiplicative noise), where (Fij(t)) = and (Fij(t)Fki(t r )) = S(t — t') (SikSji — Sudjk)- The 
latter property guarantees that Cn(t) is delta-correlated in time, Cr{€) ~ (2X OCl /kB(3 2 )S(t). Here ks is 
Boltzmann's constant, (3 — l/k^T the inverse temperature, and Aoo the heat conductivity in mean field 
approximation (p -> oo), as calculated in 0. The cross-correlations, Cdr and Crd, vanish being linear in 

r n ■ 

The heat conducting DPD solid is in fact a diffusive system with static disorder (see Ref. [12]) with heat 
conductivity A = ksfl 2 / °° dtCcjit). We apply the mode coupling theory developed for such systems in 
Rcf. [12] to calculate the LTT of Cq(t). To do so the local concentration, c(r, t) in the fluctuating diffusion 
equation of Ref. [12] needs to be replaced by the energy density e(r,i), to obtain 

d t e(r, t) = V • A(r, X) ■ V[e(r, t)/^(r, X)}, (7) 

where ip(r,X) = Cn(r) and C is the specific heat per DPD particle. The tensor A Q/3 (r,X) is the spatially 
fluctuating heat conductivity tensor with (A Q/3 (r, X)) — \5 a p, and the heat diffusivity is Dt = A/(V>) = 
X/Cn. Following the derivation of Ref. [12] one finds for the LTT, 

C Q (t) ~ C D (t) ~ -T m k B T 2 CA/(4irD T t) 1+d / 2 , (8) 

where A is the mean square fluctuation in the (k = 0) Fourier components of 5D a @(r, X) = <5A°^(r, X)/[Cn], 
i.e. 

A = (l/dV){SDf(r,X)SD^(r,X)) (9) 

with implied summation convention for repeated indices. The LTT ~ Q £ q d - m random solid 

has the same structure as the VACF in the Lorentz gas, and both LTT's vanish when the fixed particles are 
filling the sites of a periodic lattice [13], because A = 0. 

In the present system one has for large densities the explicit expression for the fluctuating heat diffusivity 
D a/3 (r,X) = (X/nV)^2 i<: j w(rij)rij iCt rij^, and the quantity A can be calculated, yielding the LTT at large 
densities, 
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Do the above results have any implications for dense classical fluids, say in the vicinity of the triple 
point? We propose the following scenario. In such dense systems the motion of the particles is quite 
restricted - somewhat comparable to a quenched system - and collisional transfer of momentum and energy 
is the dominant transport mechanism. Consider for example the microscopic stress tensor. If the rate 
of energy transfer between the particle pair (ij) in Qd is replaced by the rate of momentum transfer, 
i.e. the inter-particle force, Fij. y = —dV(rij)/drij. y , then Q x becomes the collisional transfer component 
S^y of the stress tensor. This suggests that in dense fluids the stress correlation function has a LTT, 
Cg(t) = V^ 1 (S' n/ (t)S' n/ (0)) ~ A c t~ 1 ~ d / 2 , similar to the velocity auto- correlation function in the Lorcntz 
gas Standard mode coupling theory for fluids only predicts that the TCF Cg(i) of the kinetic stresses, 
S* = J2i mvi X Viy, has a LTT ~ Akt~ d / 2 , @,IHEI]> and that no such tail is present in C c s (t). 

Standard MD simulations of TCF's in hard sphere fluids have not been able to establish a LTT ~ t~ d l 2 
in C|, as is consistent with theory. However, the non-equilibrium MD simulations jloj for small Lennard- 
Jones systems under constant shear rate 7 seem to suggest an (intermediate ?) LTT of Cg(t) ~ A c t~ d / 2 
with A c one or two order of magnitude larger than the value Ak, predicted by the theory for (7|. Here 
determination of the LTT involves two non-uniform limits (7 — ► 0, t — > 00). A large intermediate tail 
~ t~ d l 2 may be a crossover phenomenon for t < t cross {f{) at a small, but non- vanishing 7, where the 
collisional transfer stress correlation in dense fluids is showing the t~ d l 2 — tail of Burnett-type correlations, 
similar to those derived for Lorentz gases [l2| . 

It would be of great interest to test the proposed scenario for LTT's in TCF's in dense classical fluids by 
performing detailed computer simulations of Cq in DPD solids and fluids, as well as by developing a quan- 
titative analysis of Cq for general densities. Promising alternatives of measuring these LTT's, particularly 
in 1-D, are offered by simulating the Helfand moments fT?| . 
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